Distinguishing quantum from classical oscillations in a driven phase qubit 
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Rabi oscillations are coherent transitions in a quantum two-level system under the influence of a 
resonant drive, with a much lower frequency dependent on the perturbation amplitude. These serve 
as one of the signatures of quantum coherent evolution in mesoscopic systems. It was shown recently 
[N. Gr0nbech- Jensen and M. Cirillo, Phys. Rev. Lett. 95, 067001 (2005)] that in phase qubits 
(current-biased Josephson junctions) this effect can be mimicked by classical oscillations arising due 
to the anharmonicity of the effective potential. Nevertheless, we find qualitative differences between 
the classical and quantum effect. First, while the quantum Rabi oscillations can be produced by the 
subharmonics of the resonant frequency uio (niultiphoton processes), the classical effect also exists 
when the system is excited at the overtones, nujio. Second, the shape of the resonance is, in the 
classical case, characteristically asymmetric; while quantum resonances are described by symmetric 
Lorentzians. Third, the anharmonicity of the potential results in the negative shift of the resonant 
frequency in the classical case, in contrast to the positive Bloch-Siegert shift in the quantum case. 
We show that in the relevant range of parameters these features allow to confidently distinguish the 
bona fide Rabi oscillations from their classical Doppelganger. 

PACS numbers: 85.25.Am, 85.25.Cp 

I. INTRODUCTION 

Superconducting phase qubits [l|, Q provide a clear demonstration of quantum coherent behaviour in macroscopic 
systems. They also have a very simple design: a phase qubit is a current-biased Josephson jimction (see Fig. [Ija)), 
and its working states |0), |1) are the two lowest metastable energy levels i?o,i in a local minimum of the washboard 
potential. The transitions between these levels are produced by applying an RF signal at a resonant frequency 
uJiQ = {El — Eo)/h = AE/h. The readout utilizes the fact that the decay of a metastable state of the system produces 
an observable reaction: a voltage spike in the junction or a flux change in a coupled dc SQUID. In the three-level 
readout scheme (Fig. [T]^b)) both |0) and |1) have neghgible decay rates. A pulse at a frequency 0^21 = {E2 — Ei)/h 
transfers the probability amplitude from the state |1) to the fast-decaying state |2). Its decay corresponds to a single- 
shot measurement of the qubit in state Alternatively, instead of an RF readout pulse one can apply a dc pulse, 
which increases the decay rate of |1). 
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FIG. 1: Phase qubit (a) and its Josephson energy (b). The metastable states |0) and |1) can be used as qubit states. 

One of the effects observed in driven phase qubits is Rabi oscillations 0, Q : coherent transitions in a quantum 
two-level system under the influence of a resonant perturbation, with a much lower frequency dependence on the 
perturbation amplitude A via SIr = \/ -\- 5"^ , where 6 = uj — wio is the detuning of the driving frequency lo from 
the resonance frequency luiq. In resonance, Or = A. Multiphoton Rabi oscillations, at loiq = nuj {n stands for an 
integer), observed in such qubits 0, [l], were also interpreted as a signature of coherence. 

The quantum coherent dynamics in phase qubits was tested by several complementary methods (see, e.g., Refs. 
0, SSB])- In particular, in Ref. Q Rabi oscillations between the ground and excited state of the phase qubit were 
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measured by applying a 25 ns pulse at ujiq followed by a measurement pulse at u>2i- The probability to find the system 
in the upper state oscillated with the amplitude of the resonant signal, as it should in case of Rabi oscillations. In 
Rcf. 0] Rabi oscillations were seen instead as the resonant pulse duration varied; a dc readout was used. 

However, it is known that the quantum mechanical behaviour of a quantum two-level system can still be similar 
to the dynamical behaviour of the classical non-linear oscillator [1, [l^, [H!, Hi [3- particular, it was recently 
pointed out [3, flU that due to the nonlinear behaviour of current-biased Josephson junctions, a similar effect can 
arise in a purely classical way. Though direct tests - a la Bell - of whether a given system is quantum or not are 
possible, even their simplified versions (e.g., Ref. (isj ) are demanding. This motivated us to further investigate the 
classical behaviour of a phase qubit and find that there is a possibility to distinguish the quantum Rabi oscillations 
from its classical double, by the shape of the resonance, by the fact that the classical effect can be also produced by 
the overtones, nwio, of the resonance frequency, and by the sign of the resonant frequency shift. (An observation that 
a non-Lorentzian shape of the resonance should exclude a classical explanation was made already in Ref. , where 
the spectroscopy of two coupled qubits was performed. Also, a symmetric versus asymmetric Stark shift in a qubit 
playing the role of a detector was proposed to distinguish the classical and quantum behaviour of a nanomechanical 
oscillator [lH.) Classical and quantum resonances, as a function of applied drives are also studied in Ref. [l7l |. 



II. MODEL 

The phase qubit [U, Q is a current-biased Josephson junction. The Josephson potential, as a function of the phase 
difference </>, 

f/j(0) ="^(cOS0 + 0/dc//c), (1) 

forms local minima, in which the quantized metastable levels serve as the working states of the qubit. Here Ic is 
the critical current of the junction, and /dc is the static bias current. A perturbation can be produced by applying 
a time-dependent bias current, Jac sin In the quantum case, the system can be reduced to a two-level model, 
described by the Hamiltonian [0] 

A MacSintjt 

H=—a. + ^==ia^+Xa.), (2) 

where az^x are Pauli matrices, C is the capacitance of the jimction, and x ~ 1/4 in the relevant range of parameters. 
One can obtain from here the Rabi oscillations (coherent oscillations of the probability to find the system in the 
upper/lower state with the frequency fi/j) when excited near the resonance or at its subharmonics; the shape of the 
resonance is a symmetric Lorcntzian, as determined by, e.g., the average energy of the system versus the driving 
frequency uj. 

Unlike the flux qubit [H, [l^ , where the interlevel distance AE is determined by the tunneling, here AE is close 
to the "plasma" frequency ojp = [2eIc/hCY^^ of small oscillations near the local minima of the potential, Eq. ([T]): 
characteristically AE/h « 0.95a;p The same frequency determines the resonance in the system in the classical 
regime. This means that, in principle, the same AC signal could cause either the Rabi oscillations or their classical 
double fil. 



III. CLASSICAL REGIME 



In the classical regime, the phase qubit can be described by the RSJJ (resistively shunted Josephson junction) model 
[20I. [2ll|. in which the equation of motion for the superconducting phase difference across the junction, characterized 
by the normal (quasiparticle) resistance R, reads 



hC(f(j) h d(j) 
Introducing the dimensionless variables. 



we obtain 



-I- Ic sin (f) = /dc + /ac sin ujt. (3) 



■ ijjpt, 7 = a;/wp, (4) 



4> + acj) + sin (j> — r] + e sin 7r, (5) 
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where 

Idc lac 

and the dot stands for the derivative with respect to r. The sohition is sought in the phase-locked Ansatz, 

0(t) = 0o + ^(r), ^<1. (7) 

We substitute Eq. ^ in Eq. ([5|) and expand sin to third order, which yields 

•• • . sin 00 ,2 cos 00 ,T • /ON 
+ ayj + ip cos (po = rj — sm (pQ -\ — ?/) H — -i/; + e sm 7T. (8) 

2 6 

Therefore sin 0o = ''^^'^ introducing 

70 = [1-^7']'/', (9) 

we obtain 

2 

Vi + a'i/' + 7^V' = esin7T+|V^ + (10) 

which describes the anharmonic driven oscillator [2^ . 

Here we briefly point out several features of the solution of Eq. pop (for more details see Chap. 5 in Ref. 

(i) The anharmonic driven oscillator, described by Eq. (|10[) . is resonantly excited at any frequency 7 ~ ^70, where 
p and q are integers. This however happens in higher approximation in the driven amplitude e. At small amplitude 
€ the most pronounced resonances appear at 7 w 70 (main resonance), 7 w 70/2 (anharmonic-typc resonance), and 
7 « 270 (parametric- type resonance). 

(ii) The amplitude b of the small driven oscillations at the main resonance is 6max = e/(7oQ^)- When this amplitude 
is not small, then the phase-locked Ansatz becomes invalid. Then the solution of Eq. ([5]) describes the escape from 
the phase-locked state, which means the appearance of a non-zero average voltage on the contact. This voltage is 

proportional to the average derivative of the phase, </>. 

(iii) The position of the resonances for small oscillations is shifted due to the anharmonicity of the potential. This, 
for the main resonance 7ics, is given by: 

7rcs - 70 = ^^^7 (11) 

3 + 2772 



487o' 

Note that the resonance shift is negative (e.g. as in [HI). 



(12) 



(iv) The shape of the resonances, as a function of the driven frequency 7, is essentially non- symmetrical. The 
asymmetry of the main resonance becomes pronounced at 

at small enough damping. 

(v) The parametric- type resonance at 7 « 270 appears when the damping is sufficiently low, i.e. at 



o < (14) 



At the relevant parameters, r] ~ 0.9, 70 — 0.6, this means the following: 



^>^^2, (15) 
70a rj 
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which is fulfilled (see (ii)) when the solution close to the main resonance corresponds to the escape from the 
phase-locked state. 

We note that for the anharmonic driven oscillator, described by Eq. pU]) . both the resonances at 7o/2 and 270 
appear due to the anharmonicity of the potential energy and are of the same order. For the driven flux qubit in the 
classical regime [l^ the equation for the phase variable 9{t) = 60 + i]:{t) can also be expanded for small oscillations 
about the value 9q] restricting ourselves here to the linear in ^ terms, the equation can be rewritten in the form: 

ip + + 7q[1 — /isin7r]?/' = esin7T. (16) 

In this case the genuine parametric resonance at 7 « 270 appears due to the term containing sin7r • ij) (2^ . This 
explains the prevailing of this resonance over the resonance at 7 w 7o/2, due to the small anharmonicity of the 
potential in Ref . [ll] . 




FIG. 2: (Color online). Rabi-type oscillations in current-biased Josephson junctions: (a) and (b) show the time-dependence of 
the phase difference (j) and of the energy H, (c) presents the time-dependence of the energy, averaged over the fast oscillations 
with period 2-k/ui. All energies are shifted by their stationary value: Ho = 1 — 1^1 — ry^ — ryarcsinry. The parameters for 
the blue, red, and green curves are: rj — 0.95, a — 10"'', e — 10"'', and 7 = 70; for the other curves in (c) and (d) only 
one parameter was different from the above, for comparison. Namely: (c) dashed black line a — 10^'*, (d) solid violet line 
e = 2 ■ 10"^, dash-dotted black line rj = 0.9. 

Now we proceed to numerically solve the equation of motion ^ for the relevant set of parameters close to the 
experimental case. We also investigate the behaviour of the energy of the system [2l| 

H = -0^ + 1 — cos — (?7 -f- e sin7T)(/), (17) 

which determines the thermally activated escape probability from the local minimum of the potential [H] in Eq. H]). 
The classical Rabi-like oscillations are displayed in Fig. [21 In Fig. [D^a) the modulated transient oscillations of the 
phase difference </> are plotted. These oscillations result in the oscillating behaviour of the energy of the system as 
shown in Fig. [DJb). Averaging over fast oscillations, we plot in Fig. ^c) with green solid curve the damped oscillations 
of the energy, analogous to the quantum Rabi oscillations [l3| • These curves are plotted for the following set of the 
parameters: 77 = 0.95, a = 10"'^, e = 10"'^, and 7 = 70. For comparison we also plotted the energy averaged over 
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the fast oscillations for different parameters, changing one of these parameters and leaving the others the same. The 
dashed black curve in Fig. ^c) is for the smaller damping, a = 10~^; the solid (violet) line and the dash-dotted 
(black) line in Fig. EJd) demonstrate the change in the frequency and the amplitude of the oscillations respectively 
for e = 2 • 10^"^ and 77 = 0.9. We notice that the effect analogous to the classical Rabi oscillations exist in a wide range 
of parameters. 

In Fig. [3] the effect of the driving current on the time-averaged energy of the system is shown for different driving 
amplitudes: Fig. [2Ia) for weaker amplitudes, close to the main resonance, to show the asymmetry and negative shift 
of the resonance; and Fig. [3]Jb) for stronger amplitudes, to show the resonances at 7o/2 and 270 (which are also 
shown closer in the insets). We note that the parametric- type resonance at 270 originates from the third-order terms 
when the solution of the equation for tp is sought by iterations [l^l ; when there are two or more terms responsible for 
this resonance, the respective resonance may become splitted, which is visible in Fig. [3{b) for the lowest curve. An 
analogous tiny splitting of the resonance was obtained for the driven flux qubit in Fig. 4 of Ref. (loj . 
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FIG. 3; (Color online). The time-averaged energy H — Ho versus reduced frequency for relatively weak (a) and strong (b) 
driving. Different values of the driving amplitude e (multiplied by 10^) are shown by the numbers next to the curves. The 
parameters are: 77 = 0.95 and a = 10"'*. In (b) the region between the vertical black lines corresponds to the escape from the 
phase-locked state. 



6 



IV. QUANTUM REGIME 

In the quantum regime, the phase qubit can be described by the Bloch equations for the density matrix components. 
In order to take into account the relaxation and dephasing processes, the corresponding rates F^eiax = ■ -^roiax 
and = hjjjp ■ are included in the Liouville equation phenomenologically (2^ . Then the evolution of the reduced 
density matrix p, taken in the form 



1 + Z X-iY 
X + iY 1-Z 



is described by the Bloch equations (2J, |2£ 



X — — CY — A^^, 

Y = -AZ + CX - A^y, 

Z = AY- Arelax (Z - Zo) . 
where A and C stand for the off-diagonal and diagonal parts of the dimensionless Hamiltonian: 

AE 



H 



hls,c sin Lut \ ^ C* 



From these equations we obtain Z{t), which defines the occupation probability of the upper level, 

P(r)=P22(r) = i(l-Z(r)). 



(18) 

(19) 
(20) 
(21) 

(22) 
(23) 



We choose the initial condition to be ^(0) — Y{0) = 0, Z{Q) = 1, which corresponds to the system being in the 
ground state; we also consider the zero-temperature limit in which the equilibrium value of Z is = 1. 




1.5 2.0 



FIG. 4; (Color online). The time-averaged probability P of the upper level to be occupied versus the driving frequency. The 
parameters used here are: 77 — 0.95, Ej/hujp — 300, Trcis.-xi/fi'^p = V^/huj-p = 3 ■ 10"'*. Numbers next to the curves stand for e 
multiplied by 10^. Upper inset: the time dependence of the probability P{t). Right inset: the shift of the principal resonance 
(at hu ^ AE), where Au ^ - AE/h. 
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When the system is driven close to resonance, w ~ Ai?, the upper level occupation probability P(r) exhibits 
Rabi oscillations. The damped Rabi oscillations are demonstrated in the upper inset in Fig. 31 which is analogous 
to the classical oscillations presented in Fig. [21 After averaging the time dependent probability, we plot it versus 
frequency in Fig. [Hfor two values of the amplitude, demonstrating the multiphoton resonances. FigurejUdemonstrates 
the following features of the multiphoton resonances in the quantum case: (a) in contrast to the classical case, the 
resonances appear only at the subharmonics, at fia> « /S.E /n] (b) the resonances have Lorentzian shapes (as opposed 
to the classical asymmetric resonances); (c) with increasing the driving amplitude the resonances shift to the higher 
frequencies - the Bloch-Sicgert shift, which has the opposite sign from its classical counterpart. The Bloch-Sicgert 
shift (the shift of the principal resonance &i fiio k, Ai?) is plotted numerically in the right inset in Fig. |4l Analogous 
shifts of the positions of the resonances were recently observed experimentally Q . 

V. CONCLUSIONS 

In conclusion, the following criteria can be proposed to distinguish classical Rabi-type oscillations in current-biased 
Josephson junctions: (1) the appearance of resonances both at the overtones of the main resonant frequency (e.g., 270) 
and its subharmonic (70/2); (2) the asymmetric shape of the resonances; (3) a negative shift of the resonant frequency 
when increasing the driving amplitude. In recent publications these features were not reported: the multiphoton 
resonances were observed 0, [sj; the resonances observed have a Lorentzian shape 0, [l6l . [26l | and have a positive 
frequency shift which means that the observed resonant excitations in the system were in the quantum regime. 
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